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Abstract 

The two- and three-dimensional transverse-field Ising models with ferromagnetic exchange inter- 
actions are analyzed by means of the real-space renormalization group method. The basic strategy 
is a generalization of a method developed for the one-dimensional case, which exploits the exact 
invariance of the model under renormalization and is known to give the exact values of the critical 
point and critical exponent u. The resulting values of the critical exponent v in two and three 
dimensions are in good agreement with those for the classical Ising model in three and four di- 
mensions. This is the first example in which a real-space renormalization group on (2 + 1)- and 
(3 + l)-dimensional Bravais lattices yields accurate estimates of the critical exponents. 
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I. INTRODUCTION 



The real-space renormalization group framework has been developed some time ago and 
is often considered a crude approximation in practice because the truncation that is neces- 
sarily involved for tractability usually leads to unreliable estimates of critical exponents 1|. 
Standard approaches on the basis of the block-spin transformation for quantum systems 2|- 
0] succeeded in describing qualitative properties but still had difficulties in quantitatively 
accurate calculations. For the one-dimensional transverse-field Ising model, which is equiva- 



lent to the two-dimensional classical Ising model, Fernandez-Pacheco [8 1 modified the block 

n 

construction of the standard block-spin transformation [2i] to preserve the high symmetry of 
the model and could reproduce the exact values of the critical point and critical exponent 
V. Although current activities in real-space renormalization group approaches to quantum 
systems are often focused on numerically accurate evaluations, by using for instance the 
density-matrix renormalization group [3] or the multiscale entanglement renormalization 



ansatz 



10|, it is important to develop analytical or quasi-analytical methods applicable to 



higher dimensions since numerical methods are not always suitable for the calculation of 
critical properties. 



The present paper reports on our successful generalization of the one-dimensional method 
of Fernandez-Pacheco j8| to two and three dimensions. Although it is not possible to yield 
exact solutions for those higher-dimensional systems corresponding to three and four di- 
mensions in the classical representation, at least the results for the critical exponent u 
are impressive, given that quantitatively reliable real-space renormalization group methods 
in three- and four-dimensional classical models have not necessarily been established, see 
e.g. and Q-Q- 



In Section [ITl we explain the application of the real-space renormalization group method 
to the one- dimensional transverse-field Ising model and calculate the transition point and 
the critical exponents analytically. The main results of the paper are presented in Sections 
mil and IIV[ Here, we explain how to extend the real-space renormalization group method 
to higher spatial dimensions. Our study is concluded in Section |Vl 
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II. REAL-SPACE RENORMALIZATION GROUP FOR THE ONE-DIMENSIONAL 
TRANSVERSE-FIELD ISING MODEL 



Let us first review the one- dimensional scheme proposed in 8| and generahze it to the 
case with a longitudinal field, 

N~l N N 

i=l i=l 1=1 

where free boundary conditions are assumed. The system is composed of = 1/2 spins, and 
crj^, are the Pauli matrices. We consider ferromagnetic interactions J > and assume 
that hx and hz are not negative without loss of generality. This model is known to be exactly 
solvable when = 0, in which case the system is self-dual 

We start by dividing the lattice into blocks of two spins as shown in Fig. [1] The Hamil- 
tonian is also split into the intra-block and the inter-block parts. 



intra,/ 



inter, (/,/+l) 



(2) 
(3) 



where spins i and i + 1 belong to block /, and spin i + 2 belongs to block / + 1. Most 
importantly, this particular block partition is suited to preserve the form of the Hamiltonian 
under the renormalization-group transformations. 
The eigenvalues of -f^intra are 



£i - 

£3 = \J{J-hzY + hl, 
The corresponding eigenvectors are 

|1) = ai,i| tt) +a-i,i| it), 



■^{J + hzY + hl, e2 = -^JiJ-hzf + 



£2 

£4 



|2) = ai,_i| ii) +a_i,-i| U), 



|3) = a_i,_i| U) - ai,-i| U), |4) = a_i,i| ft) - Oi.il it), 



(4) 
(5) 

(6) 
(7) 



7+1 



00 00 



FIG. 1. Construction of block spins in one dimension. 
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where 

and {| tt)) I ti)) I it)> I ii)} is the orthonormal basis in the cr^-basis, i.e. t) = | t); 
^1i) = -|i). 

We next keep the two lowest lying energy eigenstates |1) and |2), and drop the others, 
|3) and |4). This procedure is expected to be effective for the study of the ground state. We 
then replace each block with a single spin representing the |1) and |2) states. To this end, 
we define the projector onto the coarse-grained system as 

N/2 

P = ®Pi, (9) 

7=1 

where Pj is the projector, 

P, = (|l)(l| + |2)(2|),. (10) 

The resulting coarse-grained Hamiltonian is PHP, and the renormalized intra-block Hamil- 
tonian is represented as 

PiHi^tre.,iPi = ^ (^1 + ^2) 1/ + ^ {si - 62) a|. (11) 
The corresponding projection of the terms in the inter-block Hamiltonian is written as 



Pj{U^al,)Pj = al (12) 

(13) 



Pi+i ® li+s) Pi+1 = ^ {all - al^i - a\i + a^i,_i) 



Pi (li ® (jf+i) Pj = (ai,ia_i,_i + ai,_ia_i,i) . (14) 
The renormalized Hamiltonian is expressed as 

N/2 N/2~l N/2 N/2 

PHP = J2 '-^h -jYI ^i^hi -h.Y. - E (15) 

7=1 1=1 1=1 1=1 

where 

J =1^ + - - a\-i) ) (16) 
K = K (ai,ia_i _i + ai _ia_i,i) , (17) 
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TABLE I. The transition point fcc,xi the exponents yx (related to h^) and Uz (related to hz) of 
the linearized renormalization group transformation for the one-dimensional transverse-field Ising 
model. The critical exponents and rj are determined from u = 1/yx and rj = d — 2yz + 2, 
where d (= 2) is the spatial dimension of the corresponding classical system. The exact solution 
is for the two-dimensional classical Ising model, which is equivalent (through a quantum-classical 
mapping to the one-dimensional transverse-field Ising model. 







Vx 


Vz 


V 


r? 


exact solution 


1 


1 


1.875 


1 


0.25 


renormalization group 


1 


1 


1.543 


1 


0.914 



hz = -- {ei - ^2) + - (a?,i - - + a\_^) + h,. (18) 

Note that our transformation preserves the form of the Hamiltonian. Other choices of 
the intra- and inter-block Hamiltonians lead to more inconvenient transformations that do 
not preserve the form of the Hamiltonian. In other words, our method does not generate 
additional coupling constants under renormalization. 

Equations (fT5l) -(fT8l) have a non-trivial fixed point at {kx,kz) = (1,0), where kx denotes 
hx/J and kz is for hz/J, which is the exact critical point of the transverse-field Ising model 
in the absence of longitudinal fields. Keeping the Hamiltonian for hz = self-dual through 
the transformations would account for this outstanding result. We can also calculate the 
eigenvalues of the linearized renormalization group transformation, and determine the criti- 
cal exponents u and 77. The results are shown in Table [H It is remarkable that the exponent 
u is exact, since real-space renormalization group calculations seldom yield exact results. 
However, the other exponent rj, which is related to the longitudinal field, is not exact. The 
longitudinal field disturbs the spectral symmetry essential to obtain exact values for the 
critical exponents. It is natural to think that the key to obtain the exact transition point 
is the existence of self-duality, while the property essential for an accurate determination 
of critical exponents may not always coincide with the self-dual character of the problem. 
The construction of block Hamiltonians described above leads to a good estimation of the 
critical exponent u and preserves the self-duality in one dimension when hz = 0. The same 
construction is effective also in the higher-dimensional transverse-field Ising model, which is 
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FIG. 2. Construction of block spins in two dimensions. This block partition preserves the form of 
the Hamiltonian under renormalization group transformations. 

not self-dual. This fact is discussed in the subsequent sections. 

III. GENERALIZATION TO TWO DIMENSIONS 

We next generalize the method to the two-dimensional transverse-field Ising model on 
the square lattice. For simplicity, we set hz = 0, and the Hamiltonian becomes 

H = -jJ2^t^]-h.J2^^^ (19) 

where spins interact with their nearest neighbors The key idea consists of performing 

renormalization group transformations which preserve the form of the Hamiltonian by a 
projective isometry that preserves the bond algebra (i.e. the algebra realized by the oper- 
ators CTiCj and af). The method to construct blocks and block Hamiltonians is especially 
crucial. Exploiting our experience in one dimension, we divide the lattice into blocks just 
as in one dimension (Fig. [2]). Furthermore, we combine the one-dimensional block method 
in horizontal and vertical directions to restore the symmetry of the lattice: If we iterate 
the renormalization in the same way as in one dimension, the system will be rescaled in 
only one direction. To renormalize the system also in the other direction, we iterate the 
renormalization in two directions; first in the horizontal direction and then in the vertical 
direction. 

In the first step of the renormalization (in the horizontal direction) we can replace each 
block with a single spin using the same procedure as in the one-dimensional case. For = 0, 
Eqs. dH) and (jH]) are 

ei = e2 = -y^J^ + hl (20) 
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and the equations corresponding to Eqs. ffT^ . ( 1T3|) and ([H]) are 



Pj {U ® Pj = a|, (22) 
® 1.+3) = -j^==a]+,, (23) 

^/ (1. ® <i) ^/ = ^=f (24) 

We find that the z-component of the spin on the right spot in a block becomes the z- 
component of the block spin, but the z-component of the spin on the left spot in a block 
becomes the ^-component of the block spin multiplied by J / \/ J'^ ^ h\. Now, let us redefine 
the coupling constants for the horizontal direction and the vertical direction as and Jy 
in order to distinguish these two quantities in this scheme. The renormalized coupling 
constants and transverse field are then written as, 

^ ( A " • 

K = , ^ ■ (27) 

In Eq. (!26|) . JyJ\l{J\ + /i^) is derived from the coupling of two spins on the left spot in 
each block, and the rest is derived from the one on the right spot in the blocks. 

Next the system is renormalized in the vertical direction in the same way as the horizontal 
direction to recover the symmetry. The coupling constants and the transverse field are now 



4 = jj + 1 , (28) 




(29) 



^ . (30) 



Jl + hi 



A 



B 



O 

o 



o 
o 



o 
o 



FIG. 3. The order of the renormahzation in the horizontal and vertical directions. With A, we 
renormalize the system in the horizontal direction and then in the vertical direction. B has the 
reverse order. 



From these equations, the following renormalization group equations are generated, 



. klK^il + kl) {(2 + klf + klkl (1 + kl)] 

kh = o ! (31) 

2{2 + klf + klkUl + kl) 



{2 + klf 



These equations can be represented as 



h = f{k.,kl), (34) 



22/(1 + 1/) 



These renormalization group equations are still asymmetric in kh and k^. To render the 
renormalization symmetric we renormalize the system in the reverse order, vertical and then 
horizontal. We define the order of the renormalization in the horizontal direction and then 
in the vertical direction as the order A, and the reverse order as B (Fig. |3]). Hence, the 
new step is represented as AB. The renormalization group transformation with order A is 
established by Eqs. (!33|) and (134|) while with order B is 



kt = f[kt,k,), (36) 
k = f {kt, kl) . (37) 
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TABLE II. The critical point k^, eigenvalues of the renormalization group transformation A, critical 
exponent i^, the slope of the dominant eigenvector 0i, and the scalar product of two eigenvectors 
01 ■ 02 from real-space renormalization group with symmetrization on the square lattice. 





A 


AB 


ABBA 


ABBABAAB 


three-dimensional classical 




1.544 


1.793 


1.828 


1.835 


3.4351 [14] 




2.383 


1.897 


1.855 


1.848 




Ai 


3.083 


9.377 


87.901 


7731.18 




A2 


1.297 


1.706 


2.912 


8.477 






0.61555 


0.61935 


0.61941 


0.61936 


0.6301(4) [15] 




2.66356 


2.59466 


2.59372 


2.59443 




slope of 01 


1.758 


1.059 


1.026 


1.002 




01 • 02 


0.644 


0.140 


0.071 


0.050 





The transformation with order AB is obtained from the substitution of kh, Eq. fl33l) . for k^ 
in Eq. ( j37l) and k^, Eq. (j34l) . for k^ in Eqs. ( l36l) and (j37l) . Although B relaxes the asymmetry 
in A, the renormalization group equations with AB is still asymmetric. The symmetrization 
procedure is therefore repeated as ABBA and ABBABAAB. 

If we regard the renormalization map in the order ABBABAAB of scaling factor 2^ as a 
single transformation, the eigenvalues of the linearized renormalization group transformation 
are 7731.18 (= {2^f-^^^^^) and 8.477 (= {2^f-^^^^^^). One of the eigenvalues is much larger 
than the other, and hence we may be justified to ignore the smaller eigenvalue. The value 
of the critical exponent v derived from the larger eigenvalue is 0.61936, which is very close 



to the reliable numerical result 0.6301 



15|- 



The effectiveness of the symmetrization scheme is clearly seen in Table [TTl With only A, 
kc^h and k^.^^ are very different, and both eigenvalues of the linearized renormalization group 
transformation are relevant. The introduction of reverse order relaxes the asymmetry of the 
fixed point. In addition, in the space of kh and ky, the slope of the dominant flow approaches 
1, and the other flow becomes perpendicular to the dominant one (Fig. S]). These changes 
suggest that the larger eigenvalue is the reliable one and we may ignore the other. The latter 
seems to be an artifact of the approximation. 

We can also study the case with a longitudinal field, 7^ 0. With the same scheme as 
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symmetrization 





FIG. 4. The renormalization group flow near the fixed point in the space of kh and k^. With the 
symmetrization procedure kcfi and kc^v become close to each other, and the slope of the dominant 
flow (the thicker line) approaches 1, and the subdominant direction becomes perpendicular. 



TABLE III. The exponents yx and t/z for the linearized renormalization group transformation and 
the critical exponents i' and rj for the square lattice, derived from symmetrization. 





A 


AB 


ABBA 


ABBABAAB 


three-dimensional classical [15j 


Vx 


1.62456 


1.61459 


1.61445 


1.61456 




Vz 


2.39774 


2.38895 








V 


0.61555 


0.61935 


0.61941 


0.61936 


0.6301(4) 


ri 


0.20452 


0.22209 






0.0364(5) 



for hz = 0, we obtain the eigenvalues of the linearized renormalization group transformation 
and critical exponents as listed in Table UTTI The exponent yz, or the critical exponent y, i s 



not very close to the numerical result of the classical Ising model in three dimensions |15| . 
a situation similar to the one-dimensional case. Moreover, this value is not improved by 
the symmetrization process. This property may give an important clue for the study of the 
relation between the symmetry of a model and the symmetry of a method. 



IV. THREE DIMENSIONS 

Next, we generalize the scheme to a cubic lattice in three dimensions. Let us again restrict 
ourselves to the system without longitudinal field for simplicity. We have to renormalize the 
three-dimensional system in three directions. The order A is defined as horizontal, then 
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FIG. 5. The order of renormalization in three dimensions. With A, we renormahze the system in 
the horizontal, then vertical and finally in the third direction. B realizes the reverse order. 



vertical and finally along the third direction. The order B is the reverse of A (Fig. [5]). We 
now define the coupling constant for the third direction as Jt in addition to Jh and for 
the horizontal and vertical directions. 

In the first step of order A, which is the renormalization in the horizontal direction, the 
parameters change as 

Jh = -F^^=, (38) 
+ (39) 
^.^..(^.l), (40) 

^ ^/jf+M' ^^^^ 

Note that the coupling constants of vertical and third directions are changed under the same 
rule. In general, when the system is renormalized in the direction a and another direction 
is represented as /3, the parameters are transformed as 

= /} ^ , (42) 
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TABLE IV. The exponents and for the hnearized renormahzation group transformation and 
the critical exponents v and for the cubic lattice, derived from symmetrization. 



AB classical mean-field 



2.0213 2.0092 2 

3.1507 3.1428 3 

y 0.49474 0.49772 i 

7? -0.30145 -0.28553 



hi 



\/ Jl + hl 

Carrying out the three steps of the scheme with these relations, we can obtain the parameters 
of the system renormalized in the three directions. 

The symmetrization with the combinations of A and B as in the two-dimensional system 
improves the result (Table II VI) . The value of v with A is 0.49474 and with AB is 0.49772. In 
the classical mean-field Ising model, which corresponds to the three-dimensional transverse- 
field Ising model, the value of v is 1/2. Our procedure yields almost the exact value of z/, and 
the symmetrization is an effective way to improve the results just as in the two-dimensional 
case. 

The result of the exponent related to a longitudinal field in three dimensions is also shown 
in Table IIVI The value of rj is not as good as the result for z/, a situation similar to the 
one- and two-dimensional cases. 



V. CONCLUSIONS 

In this paper, we proposed a real-space renormahzation group procedure for the transverse- 
field Ising model in finite dimensions. The procedure is based on the block-spin transfor- 
mation, and the preservation of the form of the Hamiltonian under the transformation is 
essential to successfully extract the critical point and the critical exponents. A remarkable 
feature of the method is that it yields exact results in one dimension js]. We have general- 
ized this idea to higher dimensions. The one- dimensional block method is also effective in 
higher dimensions, and we have combined the method in horizontal, vertical and the third 
directions to restore the symmetry of the lattice. Our results demonstrate the utility of the 

12 



block Hamiltonian we have adopted. Although the results fall short of the exact solutions 
as in the one-dimensional case, they still represent important steps because the real-space 
renormalization group rarely gives accurate estimates of critical exponents in three and four 
dimensions. 

There are several points to be clarified concerning the present method. We are particularly 
interested in establishing the reasons why the one-dimensional case yields exact results since 
the answer may give an important hint on how to improve the higher-dimensional cases. The 
exact result in one dimension strongly suggests that one of the key points is the fact that the 
one-dimensional transverse-field Ising model is self-dual. Since the critical exponent v has 
been estimated to good accuracy in two and three dimensions where there is no self-duality, 
additional factors should have contributed such as the preservation of the bond algebra 
embedded in our construction of the block Hamiltonian. It is necessary to clarify what 
has been the essential ingredient for further developments. It is also a topic of interest to 
generalize our technique to other systems including the Potts model and disordered systems. 
The latter is important due to the lack of reliable analytical approaches to three-dimensional 
systems. 
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